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ABSTRACT 



o 

, Context. Since 1997, BL Lacertae has undergone a phase of high optical activity, with the occurrence of several prominent outbursts. Starting from 
Oh 1999, the Whole Earth Blazar Telescope (WEBT) consortium has organized various multifrequency campaigns on this blazar, collecting tens of 
J. thousands of data points. One of the main issues in the study of this huge dataset has been the search for correlations between the optical and radio 
£^ , flux variations, and for possible periodicities in the light curves. The analysis of the data assembled during the first four campaigns (comprising 

■ also archival data to cover the period 1968-2003) revealed a fair optical-radio correlation in 1994-2003, with a delay of the hard radio events of 
^ ~ 100 days. Moreover, various statistical methods suggested the existence of a radio periodicity of ~ 8 years. 

i i ' Aims. In 2004 the WEBT started a new campaign to extend the dataset to the most recent observing seasons, in order to possibly confirm and 
better understand the previous results. 

_ , Methods. In this campaign we have collected and assembled about 1 1000 new optical observations from twenty telescopes, plus near-IR and radio 

■ data at various frequencies. Here, we perform a correlation analysis on the long-term W-band and radio light curves. 

\Q Results. In general, we confirm the ~ 100-day delay of the hard radio events with respect to the optical ones, even if longer (~ 200-300 days) 
y—( _ time lags are also found in particular periods. The radio quasi-periodicity is confirmed too, but the "period" seems to progressively lengthen from 
\^ . 7.4 to 9.3 years in the last three cycles. The optical and radio behaviour in the last forty years suggests a scenario where geometric effects play a 
y—( •■ major role. In particular, the alternation of enhanced and suppressed optical activity (accompanied by hard and soft radio events, respectively) can 
• ' be explained in terms of an emitting plasma flowing along a rotating helical path in a curved jet. 

CJ") Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: BL Lacertae - galaxies: jets - 
Q\ ' galaxies: quasars: general 

o : 
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^ 1. Introduction sensitive to the long-term brightness trend, but to strictly follow 



i ., -■ the short-term flux variations, becoming bluer when brighter. 

BL Lacertae (z = 0.0688 ± 0.0002; | Miller & Hawley ll 19771) is The authors suggested that the essentially achromatic modula- 

CG , the prototype of a class of active galactic nuclei (AGNs), the tion of the flux bafje leyel on long time scales . g due tQ a variation 

BL Lac objects, which are well known for their pronounced of the relativistic Doppler beaming factor, and that this variation 

variability at all wavelengths, from the radio to the r -ray band is nkely due tQ a change of the yiewing angle In contmsU the 

(see e.g. | Villata et al. || 2002| for an extended description of this strongly chromatic fast varia bility may be the result of intrinsic 

object). It has been one of the favourite targets of the Whole energetic proce sses. The data of this WEBT campaign on BL 

Earth Blazar Telescope (WEBTB which has organized several Lac were compared to the X- ray data taken in the same period 

multifrequency campaigns on this source. Two short campaigns by BeppoSAX and RXTE by Bott cher et al.l (120031) . who found 

took place in 1999, in conjunction wit h observations of t he X- some bint of correlation 
ray satellites ASCA and BeppoSAX dRavasio et alj,l2002l) . The 

campaign carried out in May 2000 - January 2001 was particu- A subsequent campaign was carried out in 2001-2003. On 
larly success ful, with periods of exceptionally dense sampling that occasion, historical optical and radio light curves were re- 
that allowed I Villata etalj (HH to follow the intranight flux constructed back to 1968, and both colour and time series anal- 
variations in detail. The optical spectrum was found to be weakly yses were performed. Colour an alysis on a longer time pe- 

riod c onfirmed the conclus ions by IVillata et al.l (120021) . and al- 

Send offprint requests to: M. Villata lowed IVillata et al.l (l2004bl) to quantify the degree of chroma- 

* The radio-to-optical data presented in this paper are stored in the tism of both the short-term and long-term variability compo- 

WEBT archive; for questions regarding their availability, please contact nents. The optical spe ctral behaviour was further analysed by 

the WEBT President Massimo Villata (villata@oato . inaf . it). iPapadakis et al.l (12007b . who interpreted the bluer-when-brighter 

1 http://www.oato. inaf . it/blazars/webt/ mild chromatism of the long-term variations in terms of Doppler 
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factor variations due to changes in the viewing angle of a curved 
and inhomogeneous jet. 

The time series analysis performed by Villat a et"af] (|2004a) 
show ed that the main radio outbursts repeat every ~ 8 years (see 
also lCiaramella et aLll2004l) . with a possible progressive stretch- 
ing of the period. Moreover, when considering the best-sampled 
time interval 1994-2003, the optical light curve was found to 
correlate with the radio h ardness ratios , with a radio time delay 
of about 100 days. Thus. PViilata e t al. (2004a) proposed a sce- 
nario where the variability mechanism propagates downstream 
in the jet, crossing less and less opaque regions. In the inner re- 
gions it produces connected optical and hard radio events, while 
in outer zones it gives rise to softer events that are apparently 
uncorrelated with the former o nes. The WEBT light curves were 
then used by iBach e t al. (2006) for a comparison with the time 
evolution of the VLBI core and jet structure on parsec and sub- 
parsec scales. The authors showed that sometimes prominent jet 
features can become as bright as, or even brighter than, the core, 
which complicates the cross-correlation analysis between the ra- 
dio light curves and those at other frequencies. However, they 
found that the radio hardness ratios derived from the WEBT light 
curves can trace the variability of the VLBI core . Their anal- 
ysis the n confirmed the optical-radio time lags of I Villata et alj 
(2004a). Jet bending was invoked to explain both the modest op- 
tical variability in the period 1981-1996 with suppressed radio- 
optical correlation, and the low-emission region in the jet around 
1 mas from the core. 

In this paper we present optical (R band) and radio data taken 
during the WEBT campaign covering the 2002-2005 periocQ. 
Their addition to the historical light curves allows us to fur- 
ther investigate both the radio quasi-periodicity and the corre- 
lation between opti cal flux and radio hardness ratio found by 
IVillataetal.l(l2004ah . 

A new campaign on BL Lacertae was organized by the 
WEBT in the 2007-2008 observing season, including three 
pointings by the XMM-Newton satellite as well as spectroscopic 
observations with the Telescopio Nazionale Galileo (TNG). The 
results of this further WEBT effort on BL Lac will be presented 
elsewhere (Raiteri et al., in preparation). 

2. Observations and results 

During the 2002-2005 WEBT campaign, we assembled about 
11000 observations in the optical UBVRI bands from 16 
observatories (20 telescopes); in order of longitude they 
are: Yunnan, Mt. Maidanak, Abastumani, Crimean, Tuorla, 
MonteBoo, Armenzano, Perugia, Torino, Sabadell, Bordeaux, 
Roque de los Muchachos (KVA, Liverpool, NOT, and TNG), 
Boltwood, Capilla Peak, Lowell (Hall and Perkins), and Kitt 
Peak (SARA). JHK data were collected from 3 observato- 
ries: Campo Imperatore, TIRGO, and Roque de los Muchachos 
(NOT). Radio observations were carried out between 1 and 
43 GHz at: SAO RAS (RATAN-600fl Metsahovfl Noto, 
Medicinal MDSCC (PARTNeR), and U MRAO. Optical an d 
near-IR data were processed as explained in lVillata et alj (|2002). 
Radio data were collected as already calibrated flux densities 

2 The radio-to-optical data presented in this paper are stored in the 
WEBT archive; for questions regarding their availability, please con- 
tact the WEBT President Massimo Villata (villata@oato . inaf . it). 
More information at http: //www. oato . inaf . it/blazars/webt/ 

3 See e.g. Korolkov & Pa riiskiil(fl97^) : lKovalev et all j!999h . 

4 See e.g. Terasranta et al. (2005) and references therein. 

5 Seee.g. lBachetal.ld2007h . 



and complemented by data from the VLA/VLBA Polarization 
Calibration Databasqj. 

The top panel of Fig.Q]shows the optical /?-band flux-density 
light curve starting from 1994 (12418 data points). It was built 
with data collected by the WEBT in five campaigns (blue dots), 
including the 2002-2005 one, as well as with literature datcQ (red 
dots). The flux densities have been corrected for the Galactic 
extinction an d the host ga l axy co ntribution has been subtracte d 
according to IVillata et all (2002, see also Villataetajj|2004b). 
Radio light curves at 22, 14.5, and 5 GHz are shown in the fol- 
lowing three panels, together with cubic spline interpolations 
through the 30-day binned data. These splines are reported in 
the last two panels, to compar e the 22 and 14.5 G Hz variations to 
the 5 GHz ones. Following IVillata et all ((2004a), we define the 
hardness ratio H22 as the ratio between the 22 and 5 GHz flux 
densities. Analogously, Hu is the ratio between the 14.5 and 5 
GHz flux densities. The events where these hardness ratios are 
greater than 1 are highlighted in yellow and called "hard" events. 
"Soft" events are characterized by higher flux densities at lower 
frequencies. 

After a period of mo derate activity (1981-1996; see Fig. 
1 in IVillata et all l2004ah . fr om 1997 the source has shown 
several prominent outbursts in the optical band. In contrast, 
the radio activity in the last 25 years does not present a cor- 
responding change of behaviour. In particular, the major ra- 
dio outburst of 1995-1996 occurred during a "quiescent" op- 
tical phase, and it was stronger at longer wavelengths. In the 
following years the reverse was observed, with outbursts ap- 
pearing enh anced at the higher radio frequencies. This be- 
haviour led IVillata et alj d2004a) to distinguish the soft ra- 
dio events, with no optical counterpart, from the hard events, 
which are correlated with the optical outbursts. Indeed, a 
cross-correlatio n analysis with the discrete correlation func- 
tion (DCF; see |Edelson & Krolikl 11988b iHufnagel & Bregman, 
1992; iPeterson et all 1 19981) method yielded a fair correlation 
(DCF ~ 0.5) between the optical fluxes and H22, with a time 
delay of the latter of about 100 days. 

One of the main purposes of the new WEBT campaign in 
2002-2005 was to check if this result would have been con- 
firmed by the data of the new period. In the top panel of Fig. Q] 
the solid line represents the average between H22 and H\n, H melin , 
shifted in time by -100 days and properly scaled. It allows us 
to recognize that most optical outbursts have their counterpart 
in the // mean curve, but with a less clear correspondence in the 
period ~ 2002-2004.5. When cross-correlating the optical light 
curve with H mem , we obtain a r esult (g reen dots in Fig. |2]l very 
similar to that of IVillata etal.1 d2004al) . with the DCF peak at 
~ 100 days slightly enhanced (~ 0.55), and a secondary peak at 
about 200 days now more pronounced. Hence, it seems that the 
hard radio events can lag behind the optical ones with two differ- 
ent time delays. To see whether a change of behaviour occurred 
at some stage, we divided the optical light curve in two peri- 
ods. The first period runs from 1994 to 2005, with the exclusion 
of the time interval where we noticed a lack of correspondence 
between optical and radio events, i.e. JD = 2452150-2453200 
(~ 2001.7-2004.5), which represents the second period. The 
blue filled circles in Fig. [2] show the DCF between the optical 
data in the first period and H me . dn . The peak at ~ 100 days ap- 
pears strongly enhanced (<: 0.7), while the second one is essen- 



6 http : //www . via . nrao . edu/astro/calib/polar/ 

7 Besides the paper s already cited inlVillata et alj d2004bh. li t eratur e 
dat a were taken fro m lHagen-Thorn et alj j2004l) . IZhang et alj j2004h . 
an d lGu et al.U 2006). 
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Fig. 1. /?-band flux-density light curve (top, mJy) and radio light curves (Jy) at 22, 14.5, and 5 GHz with their cubic spline interpola- 
tions underplotted; grey horizontal and vertical lines have been drawn to guide the eye through the flux variations; the fifth and sixth 
panels show the 22 GHz (blue line) and the 14.5 GHz (green line) splines compared to the 5 GHz (red) one; the regions where the 
higher-frequency spline exceeds the lower-frequency one have been highlighted in yellow. In the top panel the solid line represents 
the average between H22 and H14, H mean , raised to 3.5, multiplied by 10, and shifted in time by -100 days, for an easier comparison 
with the optical light curve. 



tially unchanged. On the other hand, the red empty circles show 
the radio delays in the second period: the main peaks are now 
at about 200 and 300 days. Thus, the hard radio outbursts seem 
to lag behind the optical ones by about 100 days in 1997-2001 
with a possible lengthening to ~ 200 days, while later the de- 



lay appears to increase in the range 200-300 days, until it comes 
back to the original value in 2004-2005. 

Another issue to be investigated is the quasi-periodicity (~ 8 
years) of the main radio outbursts, particularly e vident in the 
longer-wavelength light curves (see IVillata et all l2004ai) . This 
result is based on the presence of three events peaking around 
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Fig. 2. Discrete correlation functions between the optical flux 
densities and the radio hardness ratio H mean (see text) for dif- 
ferent optical time intervals: total 1994-2005 light curve (green 
dots), total period excluding the interval JD = 2452150- 
2453200 (blue filled circles), and this interval alone (red empty 
circles). The grey vertical bars indicate the typical time delays 
of 100, 200, and 300 days. 



1980, 1988, and 1996. Moreover. IVillata et alj d2004al) also no- 
ticed a possible stretching of the period: the separation between 
the first two outbursts is noticeably less than 8 years, and the fol- 
lowing time interval is slightly longer than 8 years. The occur- 
rence of the radio outburst observed in 2005 seems to confirm 
an increasing time separation between the events, the last time 
interval being of at least 9 years (see Fig.Q]). 

To quantify the effect, we auto-correlated the 8 GHz light 
curve (the best sampled one) by means of the DCF, taking the 
main radio outbursts two by two. In Fig.[3]we show the results for 
the time intervals JD = 2443000-2449000 (~ 1976.6-1993.0; 
blue filled circles), JD = 2446000-2452000 (~ 1984.8-2001.2; 
green empty circles), and JD = 2449000-2453700 (~ 1993.0- 
2005.9; red filled circles) . By considering t he centroid of the 
peak distribution (see e.g. lRaiteri et al.U2003l) . the first two out- 
bursts are found to be separated by 2701 days (~ 7.4 years), 
the second pair by 2949 days (~ 8.1 years), and the last events 
present a time lag of 3400 days (~ 9.3 years). 

3. Discussion 

In this paper we have analysed new (unpublished) as well as old 
multifrequency data collected by the WEBT consortium during 
various campaigns on BL Lacertae to address two main issues: 
i) the correlation and time delay between the optical and radio 
flux variations and ii) the possible quasi-periodicity of the main 
radio outbursts. 

We have found that the hard radio events lag behind the op- 
tical outbursts with time scales varying from about 100 to 300 
days. In particular, the longest delays occurred in the period 
2002-2004.5. A possible explanation for this variable radio de- 
lay can be found in an inhomogeneous jet scenario, where the 
mechanism producing the flux enhancement propagates down- 
stream, triggering optical and then radio outbursts as it crosses 
the corresponding emitting regions. If the jet changes its orien- 
tation with respect to the line of sight, relativistic effects pro- 
duce a variation in the observed time scales. In particular, if the 
jet portion emitting the optical-to-radio radiation has a smaller 
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Fig. 3. Auto-correlation functions on subsequent pairs of ra- 
dio outbursts at 8 GHz. Blue filled circles refer to the period 
JD = 2443000-2449000, green empty circles to JD = 2446000- 
2452000, and red filled circles to JD = 2449000-2454000. 



(larger) viewing angle, this yields a shorter (longer) delay of the 
radio events with respect to the optical ones. However, the op- 
tical light curve in the period 2002-2004.5 does not show any 
peculiar behaviour that may suggest a significant change of the 
bulk orientation of the optical emitting region. Hence, the change 
in orientation must have occurred outward in the jet, toward the 
radio region, i.e. the jet changed its curvature. Indeed, by looking 
at the radio light curves in Fig.[T| we can see that the minimum 
radio level is achieved just around 2003, and the following hard 
radio outburst of the same year is the most delayed one with 
respect to its optical counterpart of late 2002. 

Auto-correlation analysis performed on the historical 8 GHz 
light curve revealed a progressive lengthening of the time lag 
between subsequent major radio outbursts: 7.4, 8.1, and 9.3 
years. If we assume that the intrinsic period is constant, the ob- 
served stretching might be due to increased light travel time. 
Even considering that the distance between the emitting region 
and the AGN core does not vary, a change in the orientation of 
the former departing from the AGN line of sight would aug- 
ment the distance to us and hence the light travel time. High- 
resolution VLB A maps of BL Lacertae at 43 GHz presented by 
iMarscher et al.l (120081) show that the resolved jet structure may 
have typical length scales of the order of a light year. These dis- 
tances may imply much longer de-projected lengths, so that our 
hypothesis of a ~ 1-year increase of the light travel time appears 
to be reasonable. 

As a general view of the BL Lac multiwavelength emission 
behaviour in the last about forty years, we can distinguish three 
different phases. In 1968-1980 the source appears to be active 
in both the optical and radio bands, and this phase culminates 
and ends with the exceptional 1980-1981 radio outburst, which 
presented a particularly hard s pectrum in its first major peak (see 
Fig. 1 in lVillataetaIll2004al) . In the second phase, 1981-1996, 
the reverse is observed: low optical levels and variability accom- 
panied by a reduced and softer radio activity, which however cul- 
minates in the large, soft radio outburst of 1995-1997. The third 
phase seems more similar to the first one, with strongly enhanced 
optical activity and hard radio events. This behaviour may be ex- 
plained by an orientation change of the emitting jet, where the 
optical and hard-radio emitting regions are better aligned with 
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Fig. 4. Sketch of a curved jet model 
with a rotating helical path. 



the line of sight during the first and third phases, the soft-radio 
part of the jet having a larger viewing angle. The opposite for the 
second phase. 

In addition to this scenario, we also observe the radio "pe- 
riodicity", where the major periodic events maintain the hard or 
soft nature of the phase they belong to. In other words, what- 
ever the periodicity mechanism is, the 1980-1981 outburst oc- 
curred when the radio region was very well aligned, especially 
in the harder part. During the second, 1987-1988 outburst, the 
radio (especially hard) region was probably more poorly aligned, 
while the very soft 1995-1997 event occurred when the lower- 
frequency emitting regions had a smaller viewing angle. The 
last, hard radio outburst is located well inside the hard last phase. 
The radio periodicity we observe is much better seen at longer 
radio wavelengths, which are less affected by the "noise" of the 
optical and hard-radio activity, so that the relative amplitude of 
the four major outbursts is probably essentially due to the corre- 
sponding emitting region orientation. 

Inside this long-term behaviour, probably of geometric ori- 
gin, one can wonder what the nature of the single events is, 
i.e. whether the optical and correlated radio eve nts shown in 
Fig. Q] a re due to intrinsic or geometric variations. Villat a~et~aT] 
(2004b) argued that the poor chromatism of the optical out- 
bursts compared to the strong chromatism of the fast flares sug- 
gests a geometric origin of the former. In this view, the mecha- 
nism producing the flux enhancement that travels downstream 
would be the Do ppler boosting arising from a rotating heli - 



_jg rrc 

cal path (see e.g. lVillata & Raiteril [19991: lOstorero et all 12004). 



The portion of the helical path having the smallest viewing an- 
gle successively crosses and lights up the different emitting re- 
gions. Alternatively, the existence of physical disturbances going 
through the jet may also explain the observed correlated events, 
but in this case the above-mentioned different chromatism does 
not have a simple explanation. 

Figure|4]helps clarify the "geometrical" model, sketching the 
part of the jet emitting optical-to-radio frequencies. The emit- 
ting plasma flows along a rotating helical path. The yellow spot 
marks the portion of the helical path having the minimum view- 
ing angle when the line of sight is represented by the yellow 
arrow. As the helix rotates, this zone shifts along the jet, thus 
originating optical and then radio outbursts when crossing the 
corresponding emitting regions. The orientation of the curved jet 



with respect to the yellow line of sight yields a Doppler enhance- 
ment of the optical activity. In contrast, when the jet moves and 
the line of sight is represented by the blue arrow, the alignment of 
the optical emitting region changes and consequently the optical 
activity is suppressed. At the same time the soft-radio emitting 
region now has a smaller viewing angle and hence the soft-radio 
activity increases. This bulk jet motion would thus explain the 
alternate epochs of optical (and hard-radio) and soft-radio activ- 
ity noticed in the historical light curves. 
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